APERY LIMITS AND SPECIAL VALUES OF ^-FUNCTIONS 

YIFAN YANG 

Abstract. We describe a general method to determine the Apery limits of 
a differential equation that have a modular-function origin. As a by-product 
of our analysis, we discover a family of identities involving the special values 
of L-functions associated with modular forms. The proof of these identities is 
independent of differential equations and Apery limits. 



1. Introduction 

In 1978, R. Apery proved that £(3) = X^^Li n ~ 3 ^ s an irrational number by 
constructing two sequences 




and then showing that b n /a n converges to £(3) fast enough to ensure irrationality 
of C(3) (see [S]). Another remarkable discovery of Apery is that a n and b n satisfy 
the recursive relation 

(n + 2) 3 u n+2 - (34n 3 + 153?i 2 + 231n + 117)u n+ i + (n + l) 3 u n = (u n = a n or b n ). 

Thus, if we set A(t) = J2^=o a nt n an d B(t) — Yl < n=o^nt n ^ then the functions A(t) 
and B(t) satisfy the differential equations 

(1) (1 - 34t + t 2 )9 3 A + {it 2 - hlt)6 2 A + (3t 2 - 27t)9A + (t 2 - 5t)A = 
and 

(1 - 34t + t 2 )9 3 B + (M 2 - hlt)Q 2 B + (3t 2 - 21t)0B + (t 2 - 5t)B = Qt, 

where 8 denotes the differential operator td/dt. 

Apery also had an analogous result for £(2) = 7r 2 /6. He showed that if {a n } and 
{b n } are sequences of rational numbers satisfying the recursive relation 

(n + 2) 2 u n+2 - (lln 2 + 33n + 25)u n+ i — (n + 1) 2 U„ = (u n = a n or b„) 

with the initial values 

o_i = 0, a = 1, b = 0, bi = 5, 
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then b n /a n converges to £(2). Again, due to the recursive relation, the generating 
functions A(t) = ^ a n t n and B(t) = b n t n satisfy 

(2) (1 - lit - t 2 )6 2 A - {lit + 2t 2 )6A - (3t + t 2 )A 

and 

(1 - lit - t 2 )6 2 B - {lit + 2t 2 )0B - (3i + t 2 )B 

respectively. 

A common feature of Apery's two examples is the existence of a differential 
equation LA{t) — with regular singularities whose local exponents at t = are 
all such that if A{t) = 1 + ^ a n t n is the unique holomorphic solution at t = 
and B{t) = t + b n t n is the unique holomorphic solution of the inhomogeneous 
differential equation LB{t) = t at t = 0, then the ratios b n /a n converge to a special 
value of the Riemann zeta function. Inspired by these two examples, Zudilin et al. 
[H HI M considered the Apery limits of a differential equation. The general setting 
is described as follows. 

Let Lf{t) = be a linear differential equation with polynomial coefficients and 
regular singularities. Assume that the local exponents of L at t = are all so 
that the monodromy around t — is maximally unipotent. In other words, the 
differential operator L takes the form 

9 k + tFi(0) + t 2 P 2 {9) + ---+t d P d {6), = td/dt, 

where Pj are polynomials of degree < k. Let A{t) = 1 + ^a„t" be the unique 
holomorphic solution at t = 0. Then the sequence {«„} satisfy a {d + l)-term 
recursive relation 

(n + d) k a n+ d + Pi(n + d - I)a„ +d _i H h P c i(n)a n = 

with initial values a-d+i = ■ ■ ■ = a_i = and ao = 1. Now for each integer j from 
1 to d — 1, we let B(t) = P + X^j+i b n t n be a solution of the inhomogeneous 
differential equation LB(t) = j k V . (Note that when j > d, there may not exist a 
solution of LB(t) — P that is holomorphic at t — 0.) The coefficients b n also satisfy 
the recursive relation 

(n + d) k b n+d + P 1 (n + d- l)b n+d -i + ■■■ + P d {n)b„ = 

with initial values bj-d+i = • •• = bj-i = and bj = 1. Then the jth Apery limit 
is defined to be the limit of b n /a n . For example, in [8], Zudilin gave a sixth order 
differential equation whose Apery limit is the Catalan number J2^=o(~ l)"/(2n+I) 2 
and a fifth order differential equation whose Apery limit is C(4)- In [9], he also 
found a third order differential equations whose first and second Apery limits give 
simultaneous approximations to log 2 and C(2)- The paper also contains an example 
of a sixth order differential equation that gives simultaneous approximation to £(2) 
and C(3)- More recently, Almkvist et al. [1] consider the Apery limits of fourth 
order differential equations of Calabi-Yau type. The numerical computation finds 
that the recognized limits are all rational combinations of values of zeta functions 
or Dirichlet series associated with odd characters modulo 3 or 4. 

In practice, one would be interested only in the cases where the coefficients a n are 
integers having nice properties, such as having a closed form in terms of binomial 
coefficients. When such differential equations have order 2 or 3, they are often 
related to modular forms. This is because in these cases the monodromy groups 
happen to be nice arithmetic subgroups of S'L(2,R). Therefore, it is possible to 
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use the theory of modular forms and modular functions to determine the Apery 
limits for these cases. In fact, this has been done by Beukers [2] earlier. However, 
we remark that in [3j the main goal is to give a modular-function interpretation of 
Apery's irrationality results, so the argument is not readily applicable to general 
situations. Thus, the main purpose of the present article is to present a method that 
works for more general differential equations. In Section [2] we will briefly review 
Beukers' argument, and then describe our general approach. In Sections [3] and |4j 
we specialize our method to the cases where the non-zero singularity closest to the 
origin corresponds to an elliptic point or a cusp of the underlying modular curve. 
The examples we work out include cases (e), (h), and (/?) of pQ. These are some 
of the examples where [1] fails to determine their Apery limits. We also give two 
examples in which the Apery limits are values of L-functions associated with cusp 
forms. 

As a by-product, our analysis in one particular example leads us to a family of 
identities involving the special values of L-functions associated with modular forms 
of weight 3. (See Lemma[5]in Sectional) A typical example is 



for the cusp form f(r) = 77(2r) 3 7y(6r) 3 = Y^Li c ra e 2,rm ' r of weight 3. The proof 
of these identities is independent of differential equations and Apery limits. We 
expect that our argument can be extended to modular forms of higher weight, but 
we do not attempt to do so since this will be too far astray from the main topic of 
the paper. 

After finishing this paper, Zudilin kindly informed of us that in an unpublished 
manuscript 7 , Zagier also has described an approach to the determination of Apery 
limits very close to ours (presumably in the same spirit as [5]), but we do not know 
to what extent the two papers overlap. 
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Professors Yui and Zudilin for their interest in the work and providing valuable 
comments. The work was done while the author was a visiting scholar at the 
Queen's University, Canada. The author would like to thank Professor Yui and the 
Queen's University for their warm hospitality. The visit was supported by Grant 
96-2918-1-009-005 of the National Science Council of Taiwan and Discovery Grant 
of Professor Yui of the National Sciences and Engineering Council of Canada. 



Recall that a result of Stiller [1] (see also [5]) states that if A(t) is a (mero- 
morphic) modular form of weight k with character x an d t(r) is a non-constant 
modular function on a subgroup T of SL(2, R) commensurable with SL(2, Z), then 
A, tA, . . . , r k A, as functions oft, are linearly independent solutions of a (k + l)-st 
order linear differential equation LA = with algebraic functions of t as coefficients. 
We assume that the differential equation LA = has polynomial coefficients. This 
is true whenever t(r) is a uniformizcr (Hauptmodul) of the modular curve X(T). 
From now on we assume that T is of genus zero so that such a uniformizer exists. 



n=l mod 12 




Acknowledgments 



2. Modular-function approach 



4 



YIFAN YANG 



We assume that t(ioo) = and the Fourier expansion of A(t) starts from 1. 
Then the differential equation satisfied by A and t has local exponents all equal 
to at t = 0. When the function i : r h t(r) is locally one-to-one and A(t) is 
holomorphic at a point tq, the function A{t) is a well-defined analytic function of t 
near to = £(to). Conversely, if either t : t h- » £(t) is not locally one-to-one at r = tq 
or .A(to) = oo, then to is a singularity of the differential equation. Thus, singular 
points to = t(ro) of the differential equation LA(t) — can occur at four types of 
points where 

(1) to corresponds to an elliptic point of the modular curve -X'(r), or 

(2) to corresponds to a cusp of X(T), or 

(3) A(t ) — oo, or 

(4) to is a branch point of the covering map X(T) — > P X (C) given by h(r) — > 
t(r), where h(r) is a Hauptmodul of T. 

Examples. 

(1) Apery's differential equation (JTJ) has a modular-function parameterization 
given by 

( ) t[T) -{ V (2rM3r)J ' Wt)»7(6t))« " 

They are modular on To(6) + <x>6- There are four singularities 0, 17 ± 12-\/2, 
and oo. The points £ = and oo correspond to the two cusps oo and 1/2. 
The points t = 17 — 12v2 and 17 + 12^/2 correspond to the elliptic points 
t = i/y/E and r = 2/5 + i/5y/6 of order 2 fixed by (g V) and (ao -12)' 
respectively. 

(2) Apery's differential equation @ is satisfied by 

CJO 

(4) t(r)=ql[(l- q n ) 5 (f), 

n=l 

3q^ n ^ q5n—3 g5n—2 3q^ n ~ ^ 



I*) a(t) = i + V -2 + ^ q - ^ 

where (7 = e 27rlT and (f ) is the Legendre symbol. The functions are modular 
on Ti(5). The singularities 0, 00, (11 — 5v5)/2, (11 + 5v5)/2 are the values 
of t(r) at the four cusps 00, 2/5, 0, and 1/2, respectively. 
We now briefly review Beukers' argument [2]. Assume that to = 0, t\, . . . ,t m 
are singularities of the differential equation LA = with |£i| < ^1 < • ■ ■ ■ Let 
A(t) = 1 + a nt n be the t-expansion of A. In general, the radius of convergence 
of the series is |£i|. If the inhomogeneous differential equation LB(t) = t has a 
holomorphic solution B(t) = t + J2^=2 bnt n near t = 0, this series B(t) in general 
also has a radius of convergence equal to Now suppose that there is a constant 
c such that the series B(t) — cA(t) has a larger radius of convergence, i.e., such that 
t\ is no longer a singularity of B(t) — cA(t). This would mean that the sequence 
{bn/a-n} converges to c. Moreover, the larger the ratio l^/iil is, the better the rate 
of convergence is. For example, in Apery's differential equation for £(3), we have 
fi = 17 - 12^2, h = 17 + 12V2 and t 2 /h = (17 + 12\/2) 2 = 1153.999 .... Using 
the theory of modular forms, Beukers [2] showed that with the choice of constant 
c = C(3)/6, the series B(t) — cA(t) no longer has a singularity at t%. Then the 
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exceptionally large ratio £2 / £1 is sufficient to imply the irrationality of £(3). We 
now explain why the constant c is £(3)/6. 

Lemma 1. Let T be a discrete subgroup of SX(2,R) commensurable with SX(2,Z). 
Let A(t) be a modular form of weight k and t(r) be a non-constant modular function 
on r such that £(ioo) = 0. Let 

L : 9 k+1 + r k {t)9 k + ■ ■ ■ + r (t) 

be the differential operator annihilating A. Assume that g(t) is a rational function 
oft. Then a solution of the inhomogeneous differential equation LB(t) = g(t) is 
given by 

(6) B=A rU...( fU q ^m k+1 ^ 






\ t J A 

where the integration is iterated k + 1 times, and q = e 2mT . 

Remark. We did not specify the starting points and the paths of integration in 
© because different choices just give different coefficients di in general solutions 
d Q A + dirA + • • • + d k r k A + B (t), where B (t) is a fixed solution of LB(t) — 
g{t). However, when the integrand is a holomorphic modular form X)^°=i c nq n that 
vanishes at ioo, we specify the solution to be 

B 

l n \ \ In \ ^— ' I n I In 

n=l 

In particular, this is what we refer to in the examples in the next two sections. 

Proof. The proof uses the standard method of variation. Here we only prove the 
case k = 1; general cases can be proved in the same way. 
For convenience, we set 

_ qdt/dq _ qdA/dq 

Then we have 

(7) 0A = A^, 0t = -±—. 

Gi ZmGi 

We will look for two functions pi(t) and P2(t) such that B = piA+ P2TA will solve 
the differential equation LB(t) = g(t). We have 8B = Adp x + rA9p 2 + p\BA + 
P20(tA). At this point, we make an additional assumption that p\ and p 2 satisfy 

AQpi + rA9p2 = 

so that 

6B = Pi 9A + P2 0(tA). 
Differentiating again, we obtain 

6 2 B = 8 Pi 8A + 6 P 26(tA) +p 1 2 A+p 2 9 2 {TA). 
Since A and tA are solutions of L, we find 

LB = 9 2 B + ndB + r 2 B = 9 Pl 9A + 9 P2 9(tA). 
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Therefore, if p± and pi satisfy 

\ ABp x + rA6p2 = 0, 



n+e(TA)9 P2 =g(t), 

then B = p\A + p 2 rA is a solution of LB(t) = g(t). Solving the linear equations 
and using the expressions in ([7]), we find 

2irirgGi 2irigGi 



It follows that 

, /"* rgGi dt A /"* 

B(t) = -2mA / JL ^~ + 2KiTA I 



gGi dt 
~A~t' 



Making the change of variable 1 1— ► q, we obtain 

r« rgG\ dq , o . , f 5 G? 



- -2™ A J ^£l^- + 2mrA J 



A q 

(Again, although i — > q is a many-to-one mapping, different choices of branches just 
give different solutions.) Finally, applying integration by parts to the first integral, 
we conclude that 

r q gG\dq\ dq 



B«) -*/*(/* 



A q ) q 

This proves the lemma for the case k = 1. General cases can be proved in the same 
way. □ 

Lemma 2. Let all the notations be given as in Lemma[l\ The function 

^ qdt/d q y +1 g(t) 

in the above lemma is a modular form of weight k + 2 onT with character \. 

Proof. The lemma is an immediate consequence of the well-known property that 
(qdt/dq)/t is a modular form of weight 2 on T with trivial character. □ 

Example. Let t(r) and A(t) be given as in ([3]), which satisfy Apery's differential 
equation |T]) for ((3). Then we have 

= \ (E 2 (t) - 2E 2 (2t) - 3£ 2 (3t) + 6£ 2 (6t)) , 

where -^2(1") = 1 — 24 Y^Li nq n /(l — q n ). By Lemma [1] a solution of the inhomo- 
geneous solution LB = t near t — is given by 

m-mrrrtiT)*** 

Jo Jo Jo q q q 



where 

f(r) 



( qdt /dq\ 3 t 



V t J (l-34t + t 2 )A 

= ^ (Ei(T) - 28£ 4 (2t) + 63£ 4 (3r) - 36£ 4 (6r)) 

and E 4 (t) = 1 + 240 J2™=i ™V7(1 - q n ) is the normalized Eisenstein scries of 
weight 4 on SL(2,Z). 
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Since the singularity t = 17 — 12\/2 corresponds to the elliptic point r = i/v6 
fixed by ( g V ) ' ^ ne asser tion that B(t) — cA(t) , as a function of t, is not singular at 
t = 17 — 12-\/2 is equivalent to that B(t) — cA(t) is invariant under the substitution 
t — * — 1/6t. Now we have A(— l/6r) = — 6t 2 j4(t). Thus, to prove the latter 
assertion, it suffices to show that the function 

Jo Jo Jo q q q 

satisfies 

™-fi = -^>-®)- 

For this purpose, Beukers invoked the standard method of Mellin transforms. 
By the Mellin inversion formula, one can write 

1 r 2 + ica /2ttA~ s 

E( - 1/6T, = ^l,„ r(s,i(s + 3, (^) *■ 

where T(s) is the Gamma function and 

L(s) = C(s)C(« ~ 3)(1 - 28 • 2~ s + 63 • 3~ s - 36 • 6 _s ) 

is the L-function associated with the modular form /(r). We then check that the 
integrand has simple poles only at s = and s = — 2. Moving the line of integration 
to the left of Re s = —2, counting residues, making a change of variable s h- > —2 — s, 
using the functional equation 

(?f) ~ S T( S )L(s) = - (^) " 4 r(4 - s)L(4 s), 

and then applying the Mellin inversion formula again, we deduce that ((8]) indeed 
holds. (See |2j Proposition 1] for more details. See also the examples in Sections [3] 
and[U) This basically summarizes Beukers' argument. 

As we have remarked earlier, the focus of [2] is to give a modular-function proof of 
Apery's irrationality results. Therefore, the differential equations considered there 
are very well-behaved in some sense. That is, the Apery limits c of them all have 
the special property that the the functions B(t) — cA{t) no longer have singularities 
at t\. However, in many cases, we do not need such a strong property in order for 
the Apery limits to exist. For example, consider the differential equation 

9 2 A + 3t(18<9 2 + 186 + 7)A + 729t 2 (9 + if A = 0. 

It has three singularities 0, —1/27, and oo. The local exponents at these points are 
{0,0}, {—2/3,-1/3}, and {1,1}, respectively. Thus, near t = —1/27, the solution 
A(t) — 1 + ^2 a n t n has a series expansion 

A(t) =ci(t + 1/27)- 2 / 3 + ca(t + 1/27)- 1 / 3 + • • • 

for some constants c\ and C2, at least one of which is nonzero. On the other hand, 
we can show that the solution B(t) — t + ^n,t n of 

6 2 B + 3i(180 2 + 186* + 7)B + 729i 2 (6> + 1) 2 B = t 

behaves asymptotically as 

A(t)(d Q +d 1 (t + l/27) 1/3 + •••) 
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near t = —27. Therefore, we have 

B{t) - d A(t) = A(t)(di(t + 1/27) 1/3 + ■••)• 

We then can apply the Tauberian theorems to conclude that 

b n , 
efi 

a n 

To determine the exact value of do, we use the theory of modular forms. The 
detailed computation is carried out in the next section. 

In general, if the singular point t\ = t(ri) closest to the origin corresponds 
to an elliptic point t± of order n, then the local exponents at t\ take the form 
a±/n < ■ ■ ■ < a,k+i/n for some integers dj, at least one of which is relatively prime 
to n. If the (fc + l)-times iterated integral E{t) in Lemma [T] is a holomorphic 
function of r at t\ , then the same argument as that in the previous paragraph will 
imply that the Apery limit is equal to E(ti), provided that the differential equation 
does not have another singularity ti of the same modulus as t\. Likewise, if the 
singular point t\ corresponds to a cusp, then the local exponents are all equal to a 
rational number a. In this case, we have 

A(t) = (i - h) a (c k log fe (t - h) + c fe _i log k -\t - h) + ■ ••) . 

Under the same condition on E(t), we will also be able to determine the Apery 
limit, although the rate of convergence, in general, is \b n /a n — do I ^ 1/ log n, which 
is extremely slow. 

In the following two sections, we will determine the Apery limits of several ex- 
amples using the above ideas. 



jl/3 • 



3. Elliptic point cases 

In this section, we will explain in more detail how to determine the Apery limits 
when the singularity t\ closest to the origin corresponds to an elliptic point t\. 
Throughout the section, we assume that L is a discrete subgroup of SL(2,M) of 
genus zero commensurable with SX(2,Z), A(t) is a (meromorphic) modular form 
of weight k with character x on T, and t(r) is a uniformizer of the modular curve 
X(T) so that the differential equation LA(t) = satisfied by A and t has polynomial 
coefficients. We also assume that A(ioc) — 1 and t(ioo) = so that the differential 
equation has local exponents all equal to at t = 0. We let B{t) = P ' + ■ ■ ■ be the 
solution of LB(t) = j k+1 P holomorphic at t = 0, where j is a positive integer less 
than the maximal degree of the coefficients of L. For a function /, with a slight 
abuse of notations, we write /(r) if we consider / as a function of r, and write f(t) 
if we consider it as a (multi-valued) function of t. 

Let h(t) be the coefficient of 9 k+1 in the differential equation LA = 0. We assume 
that the integrand 

fqdt/dq\ k+1 P 

f{T) = {^) aW) 

inside (|6|) is a holomorphic modular form so that its L-function converges in some 
half-plane. Note that, by the assumption that t(ioo) = and A(ioo) — 1, the 
constant term of the Fourier expansion X^^Li c n<?" of /(r) is 0. According to the 
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discussion in the previous section, the Apery limit is equal to the value of the 
(k + l)-times iterated integral 

at r\. Now to evaluate J3(ti), we express E(r) using the Mellin inversion formula 

e~ x = — / T(s)x- S ds 

Z'Kl 7c— ioo 

which holds for all c > and all complex numbers x with Hex > 0, and then try 
to deduce information about E(t\) by complex analytic argument. For example, 
if /(t) is an eigenfunction with eigenvalue e of the Atkin-Lehner involution ojjy = 
( n ~o ) ^ or some positive integer N, then by an argument similar to that in the 
proof of (JSJ) , we can show that 

E(-1/Nt) = e- l {VNT)- k E{r) + (residues). 



Then we set r = i/^/N to get the value of E(i/y/W) (provided that i k e =/= 1). 

When the elliptic element that fixes T\ is not ( ^ "q 1 ) , the situation is more 
complicated. Assume that /(r) is modular on Tq(N). If we apply the Mellin 
inversion formula in a straightforward manner and write 

E (. T ) = 2^J T{s)L{ S + k + lJ) ( — J ds, 
the complex analytic argument will yield 



e \ ^ 



n=l 

where c* denote the Fourier coefficients of {-/Nr)~ k ~ 2 f(— 1/JVr) (which is a mod- 
ular form on ro(-/V) since lon normalizes Tq(N)). At this point, it is not clear how 
one should proceed to obtain information about E{j\). In such a situation, we need 
to apply the Mellin inversion formula in a different way. The key observation is the 
following simple fact. 

Lemma 3. For a matrix (™ ^) € SX(2, R) with c > and for a pair of {unctions 
g,g* : H — > C, the relation 

■<(^)=' (CT + d)V(T) 

holds for some constant e and some integer k if and only if 



S (- + -)=eH-V (---- 
c c/ V ct c 



holds. 

Proof. We have 

T + a_/a b\ (l_ _ tT 
c c \c (!/ \ ct c, 
This proves the lemma. □ 
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Now if g(r) = c ™9" an d 9*{t) — J2^=o c n1 n m * ne a bove lemma are 

holomorphic modular forms, then the lemma yields a functional equation between 
the pair of functions 

oo OO ^ 

L l s ) = y ^ e 2*ina/c L * / x = y £n e -2™nd/c 

n—l n—l 

(See Lemma 2] below.) Using the functional equation, we deduce that 
E ( — - - - | = 4r-E (- + -)+ (residues). 

\ CT C J T K \C C / 

We then make a suitable choice of r such that — 1/cr — d/ c = t / c + a/ c = t\ and 
an evaluation of E(t±) follows. 

We now prove the functional equation between L(s) and L*(s) in the following 
lemma. The lemma may have appeared somewhere in literature. However, failing 
to locate such a reference, we give a complete proof here. Notice that when g = g* 
is a modular form of weight k on r (-ZV) that is an eigenfunction of the Atkin-Lchncr 
involution wjy with eigenvalue e, the lemma gives the familiar functional equation 

l^=) S T(s)L(s,g)=(^y T(k-s)L(k-s,g) 



for the L-function L(s,g) of g. 

Lemma 4. Let T and T* be two discrete subgroups of SL(2,WL) commensurable 
with SL(2, Z) such that the cusp oo has width 1. Assume that g(r) = ^2^L c n q n 
and g*(r) — J2n°=o c n<7™ are holomorphic modular forms of weight k onT and T* , 
respectively. Suppose that 7 = (" ^) with c > is a matrix in SX(2, K) such that 

(9) g(^)=e(cT + d) k g*(T) 



^ct + d t 

for some constant e. Then the two functions 



00 

n—l n—l 

satisfy the functional equation 

(10) (—) r(s)L(s)=i k e(—) T(k-s)L*(h-s). 



Proof. For the ease of representation, here we only prove the case where the constant 
terms Co and Cq are 0. The case Co, cjj ^ can be proved by a minor modification. 
First of all, setting r = (— dr + b)/(cr — a) in ©, we obtain 



/_dr + fc\ =c -i ( _ 1)V _ fl)fc (T) 
\ ct — a J 



which by Lemma [3] is equivalent to 

(ID 9*( T --±)=e-\-*g(-^ + a - 

\C C J \ CT C 

We now consider the integral 

poo / ■ \ 

dy. 



c c 
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We have 



DC: 



-i g ( l 2 + °L ) d y = Y c n e 2ma / c / y^ e -^ny/c dy 



\ ' n— 1 ^ ^ 

Now we break the integral into two parts as usual, one from to 1 and the other 
from 1 to oo. For the integral from to 1, we make a change of variable y i — s> 1/y 
and obtain 

»i 



o V c c / Ji \ c y c J y 



Substituting the relation (JTTJ) with r = iy into the last expression, we see that 

1 y^ 9 (^ + a -) dy = i\ f>-Y (^- d -)dy 







c c 



and 



/O \ — s /"CKD / ■ \ /"OO / ■ J 



- 



V c / Jl \ cc / Jl \cc 

By the same token, we can also show that 



c c 

• « ^ / V - : //(^ • ~) '///• 

Upon the substitution s i— » fc — s into the last expression, we immediately get the 
claimed functional equation. This completes the proof. □ 

We now give two examples to illustrate our method. 

Example 1. Consider the differential equation 

(12) 9 2 A + 3<(18<9 2 + 186» + 7)A + 729i 2 (6> + 1) 2 ,4 = 0. 

This is case (h) that p] fails to identify the Apery limit. According to PQ, Arne 
Mcurman conjectures that the limit is 27r 2 /81 — L(2,X3)/2, where L(s,X3) is the 
Dirichlet L-function associates with X3 — ( 3) ■ Meurman also suggests the following 
Ramanujan-like formula 



/n\ (_ e -T/v^)n ^ 2 7T 2 1 
V3 ) n 2 (l - (-e-Vv^)") ~~ ~8T ~ 2 ^ ,X3 > 

Here we prove that these conjectures are indeed correct. 

We first note that (fT2|) is the differential equation satisfied by 



E 

n=l 



71=1 



77(r) 12 ' v ; l + 27*(r) 

where A(r) is a modular form of weight 1 on To (3) with character 

'a b\ ( d s 



(13) rf ;-V3 
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The singular points 0, oo, and —1/27 correspond to the cusps oo, 0, and the elliptic 
point 

3 + 7^3 



to = 



(> 



respectively. Using Lemma [T] we find a solution of 

6 2 B + 3t(188 2 + 18(9 + 7)B + 729t 2 ( 

is given by 



B(t) = A(t) 



qdt I dq 



dq dq 



l) z B = t 



= A(t) / 
Jo Jo 



t](3t) 9 dq dq 
v(t) 3 q q ' 



o Jo v v j + 27i) 2 q q 

The local exponents of the differential equation at —1/27 are —1/3 and —2/3. Thus, 
near t = —1/27, we have the series expansion 

A(t) = ci(t + l/27)- 2/3 + c 2 (i + l/27r 1/3 + • ■ ■ 

for some constants Ci with at least one of c\ and ci being nonzero. Now write 

Jo Jo v(r) 3 q i 

Regardless of what the constants c\ and c 2 are, we have 

B(t) - E( T0 )A(t) = A(t)(E(t) - E(r )) = A(t)(Q + d x (t + 1/27) 1 / 3 + •••)• 

At this point, the Tauberian theorems (using, for example, [3j Corollary 1.7.3]) 
already assert that if A{t) = 1 + J2 n a n t n and B{t) = t+ J2 n bnt n , then we have 



E„=i & » 



Z^—i a 



E(ro) 



1 Wfi 



< N 



-1/3 



as N — > oo. To get the stronger statement b n /a n — ► E(t ), we consider the se- 
quences a' n — a n — a n _i and b' n = b n — 6 n _i, whose generating functions are 
(1 — t)A(t) and (1 — t)B(t), respectively. Since these two functions have similar as- 
ymptotic behaviors as A(t) and B(t) near t = —1/27, the same Tauberian theorem 
implies that 



-~E(t ) 



In particular, the Apery limit is E(tq). It remains to determine the value of E(tq) = 
B((3 + V=3)/6). 

The integrand /(r) = t](3t) 9 /t](t) 3 in the definition of E(t) is in fact the Eisen- 
stein series of weight 3 associated with the cusp with the same character as ([TJ 
It can be written as 

« n 2 (g"-g 2 ") _ ~ ^ ~ //, 

/ / ^ — ff3n / / 

n=l H n=l k=l 

It follows that 

n=lfe=l v 7 

Our strategy is to use Lemma 0] to show that 

E(-1/3t + 1/3) = -£(r/3 + 2/3) + r(r) 

T 
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for some function r(r) involving r and special values of Dirichlct L- functions. Then 
we set t = ( — 1 + y/— 3)/2. With this choice of r, we have 

1 1 _ r 2 _ 3 + _ 
~37 + 3 ~~ 3 + 3 ~~ 6 ~~ T °' 

From this we can obtain the values of E{t$). 

Write the Fourier coefficients of q n in /(r) as c„. We apply Lemma 0] with 
9 = 9*= f, 

f a b\ (2 -1 N 



1 c d/ 13 -1 



4™/3 °° r p 27rm/3 



e = ,g(x) = — 1) and define 

n— 1 n— 1 

By the Mellin inversion formula, we have 

00 „ p2-7rm/3 

(-71 



E(-l/3r + l/3) = X;- 



n 

(14) 



e 

n=l 



— 2-7rm/3r 



27T« J 3 /2-ioo V 3t / 



We then move the line of integration to Re s = —5/2 and make a change of variable 
s i ► —1 — s. (By the expression of L(s) in terms of £(s) and L(s, X3) given later in 
(HHJ), we see that this is justified.) We obtain 



i r< 2+100 /2ttA s+ 

E(-1/3t+ 1/3) = (residues) + — / T(-l - s)L*(l - s) — ds. 

27TI J 3/ 2_ ioo V 3r / 

Now Lemma H] implies that 



r(5)L*(a) = -i — r(3-s)L(3- S ). 



It follows that 



r(l-s)£*(l-s) Y2tt\ 1 2s r(s + 2)X(s + 2) 



r(-i - s)l*(i - s) = -i — ^ — ^ — - = -t 



s(s + l) V 3 / s(s+l) 



/0_\ -l-2s 

= -<(yJ r( s )L(.s + 2), 

and we have 

, ,3/2+ioo /27Tt\ _S 

E(-1/3t + 1/3) = (residues) + / L(s)L(s + 2) ds. 

27T«T Jz/2-ioo V 3 * / 

By the Mellin inversion formula again, the integral in the last expression is precisely 
E(t/3 + 2/3). That is, we have 

(15) E(-1/3t + 1/3) = (residues) + -E(t/3 + 2/3). 

r 

It remains to compute the residues. 
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We have 



n 2 | - 



3 / n s k s 



n=l fc=l 

Partitioning the double sum according to the residue classes modulo 3 and simpli- 
fying, we obtain 

(16) L*(s) = ~t(s - 2)(1 - 3 3 ^)i( S , X3) + ^L(s ~ 2, X3 )C(«)(1 - 3" s ). 

Using the fact that £(s) has zeros at negative even integers and L(s,Xa) has zeros 
at negative odd integers, we see that the integrand 

r ( .)i-(. + »(£)"" -!■(.)(£ 

x (™C(*)(1 - S 1 " 8 )!^ + 2, X3 ) + XsK(s + 2)(1 - 3- 2 - 5 )^) 

in (fT4|) has poles only at s = and s = — 1. We find that the residue at s = is 

1 /3 ' 

(17) Res s=0 = - gCWCl - 3)L(2, X3 ) + Xs)C(2)(l - 1/9). 

With the evaluations C(0) = -1/2, L(0,x 3 ) = 1/3, and C(2) = tt 2 /6, we simplify 
the expression into 

Also, the residue at s = — 1 is 

- f-|cC-l)(l - 9)i(l >X3 ) + ^£'(-l,x 3 )(l - 1/3)) §£. 
Using the functional equation 

(s) r (— J^Ks) 1,(1 

for L(s,x 3 ), on e may deduce that 

£'(-i,x 3 ) = ^(2,x). 

Together with the evaluations £(— 1) = —1/12 and L(l,x 3 ) = 7i-\/3/9, we find that 
the residue at s = — 1 can be simplified as 

„ 27riv3 1 

(18) Res s= _ 1 -^- + -L(2,x). 

In summary, the above computation (fTS|) . HZ]), (|18p shows that 

25(-l/3r + 1/3) = -^£(2, x) + ~ + 1 27T ^f + -E(r/3 + 2/3). 

2r r 81 t 

Setting t = (-1 + V^3)/2, we find 

We summarize the above computation as follows. 
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Theorem 1. Let {a n } and {b n } be the sequences satisfying the recursive relation 

(n + 2) 2 u n+2 + (54n 2 + 162n + 129)n n+ i + (n + lfu n = (u n = a n or b n ) 
with the initial values do = 1, di = —21, &o = 0, and b\ = 1. Then we have 

« n- 1 / 3 



bn 2^1 
+ T. L 2 :X3j 

a n 81 2 



as n —> oo . Moreover, the evaluation 



E 

n=l 



V3 J n 2 (l - (_ e -*/V3)n) _ "8T ~ 2 ( ' X3 ' ) 



holds. 



Example 2. In this example, we will construct a differential equation whose Apery 
limits give a good rational approximation to the value L(2, /) at 2 of the L-function 
associated with the cusp form /(r) = n(r) 3 n(7T) 3 . 

The Hauptmodul of the genus-zero group To (7) + lu-j is given by ^(t) 4 /?7(7t) 4 + 
49t7(7t)7?7(t) 4 . We set 



<(r) = + 14 + 45 w J = 9 " 1092 + 4V " 1Mq4 + 

The value of t(r) at the unique cusp oo is obviously 0, and at the elliptic points i/V7, 
(5 + v 7 — 3)/14, and (7 + y 7 — 7)/14, it takes the values 1/28, 1, and oo, respectively. 

According to Lemma Q] and the discussion at the beginning of this section, we 
can basically choose any rational function git) such that <?(0) = and its numerator 
has a degree smaller than its denominator and then set 



A(r) 



f gdt/dq \ 2 g(t) 

\ t J f( T y 



Then the differential equation satisfied by t and A will give us a rational approxi- 
mation to L(2, /). 

Here we choose g(t) = t/(l — 29t + 28t 2 ). With this choice, we find 



n=l 

and the differential equation satisfied by t and A is 

(1 - <) 2 (1 - 28t)8 2 A - Ut(l - t) 2 9A - (2t + At 2 )A = 0. 

The solution B{t) of the inhomogeneous differential equation 

(1 - t) 2 (l - 28t)6 2 B - 14i(l - t) 2 9B - (2t + At 2 )B = t(l - t) 

has a ^-expansion B(t) = t + 45i 2 /4 + ■ • ■ . Furthermore, we can show that B{t) — 
L(2, f)A(t) has no singularity at t = 1/28. Therefore, we have the following good 
rational approximation to L(2, /). 

Theorem 2. Let /(r) = ?7(t) 3 ?7(7t) 3 be the cusp form of weight 3 with character 
x(cd) ~ (f) on ^o(7) andL(s,f) be its associated L-function. Let {a n } and {b n } 
be sequences satisfying 

(n+3)V l+3 -(30n 2 +134n+150)u„+2+(57n 2 + 142n-|-81)u n+ i-(14n+28n 2 )u, i = 
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with the initial values 

a a = I, ai = 2, a 2 = 24 , 6 = 0, b 1 = 1, & 2 = 45/4. 
Then we have 

h 

< 28"". 



--L(2,f) 

(In. 



4. Cusp cases 

Let t(r) and A(t) be given and assume that the integrand in ([6]) is a holomorphic 
modular form that vanishes at ioo as before. In this section we consider the cases 
where the singularity ti closest to the origin corresponds to a cusp a of the modular 
curve X(T). Pick an element in a = ("J) € SX(2,R) (preferably lying in the 
normalizer of T in SL(2, M)) such that ctoo = a. Note that the asymptotic behavior 
of .A(t) near a is determined by that of A(ar) in a neighborhood of oo. To be more 
precise, we have A(<jt) = (cr + d) k A*(r) for some modular form A*(r) of weight 
k on the group <7 _1 Lcr. Also, t*(r) = t(ar) — t\ is an algebraic function of t that 
vanishes at oo. Then for some positive integer r depending on the order of A*(t) at 
oo, A* (t) t has a power series expansion in f* . The term (cr+d) k gives a logarithmic 
factor dk log fe t* H h d . Therefore 

A (<tt) = (d k \og k t* H h d )(t*) a x (a power series in t*) 

for some rational number a. (The number a in fact is equal to the local exponent 
of the differential equation at t\.) This gives us the asymptotic behavior of A(t) 
near a. 

By the same token, to determine the asymptotic behavior of B(t) near t — t\, 
we need to consider B(ar) in a neighborhood of oo. Let /(r) be the modular 
form of weight k + 2 inside the integral in ©. We apply Lemma 2] to /(r) and 
/*(t) = (cr + d)~ k ~ 2 /(t), which is a modular form on ct _1 Lct. We then use the 
functional equation to find the asymptotics of B(t). The actual computation is 
similar to that for the elliptic point cases. We now give some examples. 

Example 3. Consider the differential equation 

9 2 A - 4t(86 2 + 86 + 3)A + 256i 2 ((9 + I) 2 A = 0, 

which is case (e) in pQ. It is conjectured that the Apery limit is L(2,x_i)/2, where 
X-i is the odd Dirichlet character modulo 4. We now show that this is indeed the 
case. 

The differential equation is the one satisfied by 

n(rf V (Ar)^ V (2r) 22 

h t j = — ztkiym — ' A y> ~ 



r;(2r) 24 ' v ; V (t) 12 t](4t) s ' 

where A(t) is a meromorphic modular form of weight 1 on To (4) with character 

(19) x(l ty=(-iTx-i(d), 

and t(r) is a modular function on To (4). By Lemma [TJ 

i ri r/ ( T )4 ??( - 4r )8 dq dq 



B{r)=A(r) / / 
Jo Jo 



r;(2r) 6 q q 

is a solution of the inhomogeneous B 2 B - 4t(88 2 + 89 + Z)B + 256t 2 (9 + l) 2 B = t. 
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The singularities 0, 1/16, oo correspond to the cusps oo, 0, and 1/2, respectively, 
and the local exponents are {0,0}, {—1/2, —1/2}, and {1, 1}. Since 



2i7 1 { T ) s T 1 (4T) 121 
we know that 

A(t) = (t - l/16)- 1/2 ( Cl log(t - 1/16) + c 2 + • • • ) 

as t approaches 1/16 for some nonzero constant c\. (The exact value of c\ is not 
important, as it suffices to know that it is not zero.) Let 

<•« ri J7 ( t )4 J7 ( 4t )8 dq dq 



E(r) 



In the following we will show that 



77(2t) 6 q q 



(20) Ei-l/Ar) = ^(2,X-l) ^ ' % ? ^ W * * 



2 v ' A ; 32t 8t7 7 ?7(2t) 6 g g 
This would imply that 

B(t) - ±L(2, X-i)^W = (* - l/16)- 1/2 (d + •••)• 
Then we deduce from the Tauberian theorems that 

1 



^( 2 >X-i) 

On 2 



< 



log n 

(Again, we shall apply the Tauberian theorems to the sequence a' n = a n — a„_i 
and b' n = b n — 6„_i with generating functions (1 — t)A(t) and (1 — t)B{t) to get the 
stronger conclusion.) We now prove (f2"U|) . 

The function t?(t) 4 77(4t) 8 /77(2t) 6 is in fact the Eisenstein series of weight 3 as- 
sociated with the cusp 1/2 with character given by (fT9|) . Its g-expansion can be 
alternatively written as 

v(t)M*t) 8 = g ( _ 1 ^_i n 2 q" 



t]{2t) 6 i + o 

v / n— 1 



2ri ' 



Then by the Mellin inversion formula 

£-l/4r=— / r#+2 — ds, 

27rl J3/2-ioo V 4r / 

where 

oo oo (_-nfc-l 

+ 2 ) = EC- 1 )"" 1 " 2 E ((2 fc-l)n)'+' = C(S)L(S + 2 ' " 2 " S) - 

Moving the path of integration to Res = —5/2 and making a change of variable 
sh — 1 — s, we obtain 

E(-l/4, ) = i £ ( 2)X -i)-^ 

i /-3/2+ioo / 2 7riV +1 
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Using the functional equations for £(s) and L(s, x-i) an d the Legendre duplication 
formula for F(s), we find the integral on the right-hand side is equal to 

/•3/2+ioo 

/ T( S )as + 2)L( S ,x-i)(l-2- s - 2 )(-mr)- s d S . 

/3/2— ioo 

By the Mcllin inversion formula again, this is equal to 



1 i 
2ttI 2t 



i [ q [ q j^f-l\n 2 q n / 2 dqdq_i f q f q ^{rf^rf dq dq 



8r Jo Jo V n J 1 - <t 1 1 8r J J Q ry(2r) 6 q q ' 
This establishes (|20[) . In summary, what we have shown is the following result. 
Theorem 3. Let {a n } and {b n } be the sequences satisfying 

(n + 2) 2 u n+2 - (32n 2 + 96n + 76)u„+i + (n + l) 2 u n = {u n = a„ or b n ) 
with the initial values do = 1, a\ = 12, bg = 0, and b\ = 1. Then 

1 



bn 1 Tin \ 
a n 2 



\ogn 



as n — > oo. 

Example 4. Let /(r) = 77(2r) 3 7y(6r) 3 be the cusp form of weight 3 on To (6) + 0J3 
with character 

*{L 9- (-"■(!)■ ^-H- 

We will find a rational approximation to 2). 
Choose 

(2D <M = ^TTmTT' *M-*gw 

r )( T ) Vi-JT) i](2T)r)(bT) 

where i(r) is modular on the smaller group To(6). The differential equation satisfied 
by t and A is 

(1 + 8i) 2 (l + 9t)9 2 A + 9i(l + 8t) 2 0A + 2t(l + I6t + 72t 2 )A = 0. 

The singularities 0, —1/9, —1/8, and 00 correspond to the cusps 00, 1/2, 1/3, and 
0, respectively. Set g(t) = t/{\ + 8t)(l + 9t) so that 

/ qdt/dq\ 2 g(t) 
{^) — = /(r) ' 

Let A(t) = 1 + T,^Li a nt n and B{t) = t - It 2 + Y^=3 M n be the solution of 

(1 + 8t) 2 (l + 9t)6 2 B + 9t{l + 8t) 2 8B + 2t(l + 16* + 72t 2 )B = t + 8t 2 

holomorphic at t — 0. We now compute the limit of b n /a n . For this purpose, we 
need to determine the behavior of 



a Jo 11 



near the cusp 1/2. Choose 

1 (Z -2 
CT =7H6 -3 

and consider E(ar). 



eT {6)+LO 3 , 
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Since /(err) = i(2y/3r — v / 3) 3 /(t), according to Lemma[H if /(r) = J2^=i °nq n , 
then the function 



L(s) = J2(-ir^ = -L(s,f) 



satisfies 



Also, we have 



n=l 

s-3 



T(s)L(s) = I -^J) r(3-a)L(3-s) 



7T 

V37 " ~ y^j 



1 \ 1 



E j ^ + — ) = — — / 3/2+100 T(s)L(s + 2) (^-) ds. 

\ 2V3r 2/ 2W3/2-00 Vv^t/ 

The integrand has simple poles at 5 = and s = — 1 with residues 

Res s=0 = L(2) = -L(2, /), Res s= _ x = -L(l) U^-) = -L(2, /). 
Using the functional equation for L(s) and the Mellin inversion formula, we find 

By Lemma El this amounts to 

E(ar) + L(2, /) = z(2 V / 3r - V3)~\E(t) + L(2, /)). 
Since A(t) satisfies 

A(ctt) = -i(2V3r - \/3)A(r), 

we have 

A(ar)(E(ar) + L(2, /)) = A(r)(B(r) + L(2, /)). 
From this we see that the function A(t)(E(t) + L(2,f)), as a function of t, is 
holomorphic at t = —1/9. Therefore, we have the following approximation of 
i(2,/). 

Theorem 4. Let /(r) = 7/(2t) 3 ?7(6t) 3 and L(s, /) 6e «is associated L-function. Let 
{a n } and {b n } be the sequences satisfying 

(n+3) 2 ii„ +3 +(25n 2 + 109n+120)w„ +2 +16(13n 2 +35n+24)M„ + i+144(2n+l) 2 M„ = 
with the initial values 

ao = 1, ai = —2, <i2 = 10, &o = 0, b\ = 1, &2 = — 7- 
Then we have 

« (8/9)™. 



-+H2J) 



Remark. Observe that the result in the above example shows that the composition 
t i ► t i ► F(t) = A(t)(E(t) + L(2, /)) is a single-valued analytic function at 
t = —1/9. This in turn implies that the function F(r) is invariant under the 
substitution t h> 71 for any 7 in the stabilizer subgroup of 1/2 in To (6). This 
observation leads us to some interesting identities for the values of L-functions. We 
will discuss these identities in the next section. 

Example 5. Consider the differential equation 

9 3 A - 8t(29 + 1){29 2 + 29 + l)A + 256t 2 ((9 + if A = 0. 
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This is equation (/?) in pQ. We now determine its Apery limit. 
The differential equation is satisfied by 

r?(r)^(4r) 16 ?7 (2r) 20 



r?(2r) 24 ' v 7 77(r) 8 77(4r) 8 ' 

They are modular on To (4). The singularities 0, 1/16, oo correspond to the cusps 
oo, 0, 1/2, respectively. The solution of the inhomogeneous differential equation 

6 3 B - 8<(2(9 + 1)(26> 2 + 29 + l)B + 256i 2 ((9 + 1) 3 B = t 

is B{t) = A(t)E(t), where 

E ^ = 5Zn F F F " 17E ^ + 1&E ^)) 

Z4U Jo Jo Jo 1 1 Q 

We have 

A(-l/4r) = -At 2 A(t). 

To determine the Apery limit, we should study how E{r) transforms under the 
substitution r ^> — I /At. 

As usual, we apply the Mellin inversion formula and write 

1 /■ 3 / 2 + i00 _ , /27ri 



27r« ./3/2-ioo V 4r 



£(-l/4r) = — / T(s)L(s + 3) ( — ) ds, 

where 



L(s) = C0)C(s - 3)(1 - 17 • 2- s + 16 ■ 2~ 2s ) 

satisfies 

n- s T(s)L(s) = 7r s - 4 r(4 - s)L(A - s). 

Moving the line of integration to Res = —7/2, counting the residues, making a 
change of variable s — > —2 — s, using the functional equation above, and invoking 
the Mellin inversion formula, we obtain 

From this we deduce the following result. 

Theorem 5. Let {a n } and {b n } be the sequences satisfying 

(n + 2) 3 u n+2 — 8(2n + 3)(2n 2 + 6n + 5)u n+ i + 256n 3 u n = (u n = a n or b n ) 
with the initial values ao = A, ai = 8, bo = 0) an d b\ = 1. Then we have 

1 



a n 16 



log n 



as n — ► oo. 
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5. Identities involving L- values 

Let the functions t(r), A(r), /(t), and E(t) be given as in Example 4. Fol- 
lowing the remark at the end of the example, we know that the function F(t) — 
A(t)(E(t) + L(2, /)) satisfies F( 7 "r) = F(t) for all integers n, where 7 = ( x 7 2 Z5) 
is a generator of the stabilizer subgroup of the cusp 1/2 in To (6). Now we have 
A(jt) = (12r - 5)A(t). Thus, the function E(t) satisfies E(jt) + L(2, /) = 
(12r - 5) _1 (£'(t) + L(2, /)). By LemmaEJ this is equivalent to 

(22) E + 1) + L(2, /) = - T (E (--L + 1) + L(2) /) 

On the other hand, we have /( 7 r) = (12r - 5) 3 /(r). By Lemma SJ letting c„ 
denote the Fourier coefficients of /(t) = X)^Li c n<7™) the pair of functions 

00 00 

n—1 n—1 

satisfy 

1) ~* T( S )L( S ) = (I) " 3 T(3 - s)L*(3 - s). 
Then, arguing as before, we have 

< 23 > £ (-t^ + b)= £ -< 2 >-; £ < 2 >-M5 + b 



Comparing ([22)1 and (|23|) . we find 

OO OO OO 

E°n c 147rw/12 _ £^_ e 107rin/12 _ _ £^ 
n—1 n—1 n—1 

Now we observe that the deduction of (|22|) and ([23]) is really independent of 
differential equations and rational approximations. Thus, we expect that identities 
of this type must exist in a general situation. 

Lemma 5. Let T — To(N) or Ti(N). Let a = p/q / 00 fie a cusp of T . Let 
/(t) = J]^Li c„g™ 6e a holomorphic modular form of weight 3 on T luit/i character 
X- Assume that /(r) vanishes at 00 and a. Then for any 7 = ( " d) wii/i c > m 
i/ie stabilizer subgroup of a inT such that x{l) — 1? ^awe 

L(2)=L*(2) = L a (2), 

where the L-functions are defined by 

00 00 00 

L( s ) = ^ ^e 2 ™"/ c , L*(s) = Y^^e- Mnd / C , L a (s) = —^ ma 

n—1 n—1 n—1 

for Re s > 3 and i/ien continued analytically to the whole complex plane. 

Proof. Let ^ be the path consisting of the vertical line from — d/c + ioo to —d/c, 
the semi-circle from —d/c to a, and the vertical line from a to a + ioo. Consider 
the integral 



(cr + d)/(r) dr. 



Since —d/c is equivalent to 00 under T, by the assumption that f(r) vanishes at 
00 and a, the integral is convergent. Therefore, it is equal to because /(t) is 
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assumed to be a holomorphic modular form. Now we make a change of variable 
t — ► 7 _1 r in the integral from —d/c to a. We have 

^ (cr + d)f(T)dr= f (c 1 - 1 r + d)f^r) dT = ^ f{r)dr. 
J -d/c Jioo {-CT + ay J too 

It follows that 

/—d/c+ico fCt+ioo 
(cT + d)f{r)dT = / (cr + d-l)f{T)dT. 
-d/c Ja 

Now recall that a stabilizer 7 of the cusp a = p/q takes the from 

1 + pqm, —p 2 m 
q 2 m 1 — pqm 1 

for some integer m. Hence, (f!M)) can be written as 

/— d/c+ioo />CK+ioo 
(cr + d)f(r) dr = I (cr — ca)f(r) dr. 
-d/c J a 

Finally, the two sides of the above identity are equal to 

ic „ , „ ic 
-L*(2), - M 2), 

respectively. This gives L*(2) — L a (2). By a similar argument, we can also show 
that L(2) — L a (2). This proves the lemma. □ 

Example 6. Consider the differential equation 

(1 - lit - t 2 )9 2 A - (lit + 2t 2 )0A - (3f + t 2 )A = 0, 

appearing in Apery's sequence for £(2), where t and A are given by ([4]) and (|5j). 
Using Lemma [TJ we find a solution for the inhomogeneous differential equation 

(1 - lit - t 2 )0 2 B - [lit + 2t 2 )6B - (3t + t 2 )B = t 

is given by B(t) = A(t)E(t), where 
(25) 

n 2 q n dq dq 



m = L I S 2 S + 2 E £ 

\n=l mod 5 n=2 mod 5 n=3 mod 5 n=4 mod 57 

In [2], Beukers mentioned that it can be shown that 



1 — q n q q 



5t + 1/V \ 5t + 1 

but did not give details. Here, we will prove this transformation formula using 
Lemma [5] 

Since A(r/(5r + 1)) = (5r + l)A(r), by Lemma[3j it suffices to prove that 

<-) K-^-i)-^-KKi + i)-^ 

Let /(r) = 53 n =i c ™9™ the modular form of weight 3 in the definition of E(t) 
in ([25]) . It satisfies /(t/(5t + 1)) = (5r + l) 3 /(r). Now, by the usual argument, 
we find 
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where L(s) = c n e 2nm ^ 5 /n s and L* — ^c n n 27rm / 5 / n 2 . Now applying Lemma 
[5] with ( 5 ° ) that fixes 0, we have 



L(2)=L*(2) = L(2,/) = E%, 

n— 1 

which is shown to be C(2)/5 in [2]. This gives the desired equality (|2l))) . 



Corollary 1. Let X)»Li c nQ n ^ e ^ e Fourier expansion of /(r) = 77(2t) 3 77(6t) 3 . 
-Let L(s, /) = X)«Li c n/n s be the L-function associated with f . Then we have 

E V^t^./), £ % = ^i(2,/). 

n=l mod 12 n=7 mod 12 

Proof. Obviously, c„ is whenever n is even. Also, using Jacobi's triple product 
identity, we have 



7 ? (2t) 3 t ? (6t) 3 = ^T(-l) n (2n + i) g ( 2 «+i) 2 /4 ^(_l)«(2n + 1) 9 3 



,3(2,1+1)74 

>.n=0 / \n=0 



From this we see that c n = whenever n = 5, 11 mod 12. Thus, letting L\ = 

E„ = lmodl2 C »/« 2 ; I* = E„^7mod 12 C "/ n2 > and L 3 = T,3\n C n/n 2 , We have 

L(2,/) = L\ + L2 + L3. Moreover, since C3„ = C3C n for all n, we have L3 = 
-L(2,/)/3. It follows that 

(27) Li+L 2 = |l(2,/). 

We now apply Lemma [5] with the matrix ( ^ ? ), which fixes the cusp 0. Lemma 
[5] gives 

00 

e 2 -/i 2 (il _ L2 ) + ^^_ z » = L( 2,/). 

n— 1 ^ ' 

Comparing the real parts, we find 

(28) L^L 2 = ^=L(2J). 



Combining lj27[) and J28J), we obtain the claimed identities. □ 

Corollary 2. Let X)^°=i c «9" ^ e ^ e Fourier expansion of /(r) = 7y(4r) 6 . Let 
L(s, /) = X)^°=i c n/n s be the L-function associated with f . Then we have 

EVE ^ = cosfL(2,f), 

n=l mod 16 n=9 mod 16 

E E % = -sin^(2,/). 

rl=5 mod 16 n=13 mod 16 

Proof. Clearly, the coefficients c„ vanishes for all n not congruent to 1 modulo 4. 
Let L i; i = 0, ... , 3 denotes J2^=4i+i mod 16 c n/ n2 - Apply Lemma[5]with the matrix 
( 16 l )' w hich fixes the cusps 0. We obtain 

e 27 "/ 16 (L - La) + e 10 ™/ 16 ^ - L 3 ) = L(2, /). 

Comparing the real parts and the imaginary parts of the two sides, we find 

cos^(L -L 2 )-sm^(L 1 -L 3 ) =L(2, /), sin ~ (L - L 2 ) + cos ~ (L x - L 3 ) = 0. 
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From this, we obtain the claimed identities. □ 
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